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Abstract 

The nonlinear stability of an oblique mode propagating in a two-dimensional com- 
pressible boundary layer is considered under the long wave-length approximation. The 
growth rate of the wave is assumed to be small so that the ideas of unsteady nonlinear 
critical layers can be applied. It is shown that the spatial / temporal evolution of the mode 
is governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity 
and density. Expressions for the linear growth rate show clearly the effects of wall heating 
and cooling, and in particular how heating destabilises the boundary layer for these long 
wavelength inviscid modes at 0(1) Mach numbers. A generalised expression for the linear 
growth rate is obtained and is shown to compare very well for a range of frequencies and 
wave-angles at moderate Mach numbers with full numerical solutions of the linear stability 
problem. The numerical solution of the nonlinear unsteady critical layer problem using a 
novel method based on Fourier decomposition and Chebychev collocation is discussed and 
some results are presented. 

1. Introduction 

One of the more important contributions to the linear stability of compressible bound- 
ary layers was the work of Lees & Lin (1946) who developed the inviscid theory and 
identified the different types of modes which exist in regions of relative subsonic and su- 
personic flow. This and subsequent modifications to the work to include viscous effects, is 
extensively described in the reviews by Reshotko (1976), and Mack (1984,1986) who also 
discusses the numerical results for the compressible linear stability problem. Recent work 
has concentrated on developing a more systematic approach, based on the the triple-deck 
type framework, and Smith(1989) discusses the lower- branch stability properties, whilst 
Gajjar & Cole (1989), Gajjar (1990), Cowley & Hall (1990), Brown & Smith (1990) consider 
the upper- branch and inviscid stability properties, with the last two papers concentrat- 
ing more on the description of the acoustic and vorticity modes at high Mach numbers. 
The nonlinear stability properties, particularly of the inviscid modes, has however received 
scant attention, although exceptions include the paper by Goldstein & Wundrow (1990) 
which extends the Cowley & Hall work to include nonlinearity and also the work of Gajjar 
& Cole (1989,1992) who consider nonlinear neutral modes in compressible boundary lay- 
ers, and Leib (1991) who has studied the evolution of modes in compressible shear layers. 
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One of the aims of this paper is to discuss the nonlinear evolution of the long wavelength 
inviscid modes in compressible boundary layers. 

It is now well known that the critical layer where the phase speed of the disturbance 
wave is equal to the local flow velocity plays an important role in the stability of many shear 
and boundary layer flows. The ideas of nonlinear critical layers were first put forward by 
Benney & Begeron (1969) and Davies (1970) and they have been widely used in numerous 
papers since then. An excellent review of the many different aspects and properties of 
linear and nonlinear critical layers is given in the cuticles by Stewartson (1981) and Maslowe 
(1986). Benney & Bergeron showed how the properties of the nonlinear critical layer could 
be used to describe nonlinear neutral waves in parallel flows. This work has recently 
been extended to compressible boundary layer flows by Gajjar & Cole (1989), (hereafter 
referred to as I), Gajjar (1990), and Gajjar & Cole (1992), who obtain various results for 
the nonlinear neutral frequencies and wavenumbers which depend on the amplitude of the 
disturbance wave. All these results are based, however, on the underlying assumption that 
the neutral modes exist and that the critical layer is of an equilibrium type. Although a 
description of the nonlinear neutral structure is interesting, and provides a solution of the 
compressible Navier- Stokes equations, of more importance is the question of whether the 
assumed nonlinear neutral structure is stable and indeed attainable. This question and the 
relevance of the equilibrium critical layer to the overall stability picture has so far not been 
resolved and needs further consideration. In shear layers the far downstream asymptotic 
structure of the non- equilibrium critical layer bears close resemblance to Benney- Bergeron 
critical layer, although there are modifications stemming from unsteady effects and outer 
diffusion layers, Goldstein & Hultgren (1988). 

Studies of the spatial evolution of waves on shear flows and in which the critical layer 
is of a non-equilibrium type have been conducted in a series of papers by Dr. Goldstein 
and his group, see for example Goldstein & Leib (1988), Goldstein & Hultgren (1988), 
Goldstein & Wundrow (1990). The nonlinearity studied in these papers is such that it 
induces a fully unsteady nonlinear critical layer problem whereas in other instances with 
weaker no nli nearity, the resulting evolution of the disturbance is described by an integro- 
differential equation of the Hickemell (1984) type, Goldstein & Leib (1990), Goldstein & 
Choi (1989), Leib (1991). 

These papers demonstrate that even though spatial equilibration can occur in some 
instances, with the properties of a quasi-equilibrium critical layer coming into play, there 
are important differences between the equilibrium and non-equilibrium type approaches, 
see also Gajjar & Smith (1985). Comparisons of the theoretical predictions and experimen- 
tal results for shear layers based on the latter approach show good agreement, Hultgren 
(1992). Mankbadi (1992ab) had also made extensive comparisons with experiments of the 
the incompressible weakly nonlinear theories for the Blasius boundary layer and his results 
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give convincing quantitative evidence that the unsteady nonlinear critical layer theories do 
indeed capture the early stages in the transition to turbulence of these flows. 

Wu (1992) has studied the weakly nonlinear instability of Stokes layers using similar 
ideas and also obtained a novel generalisation of the Hickernell integro- differential equation 
by incorporating a weak spanwise dependence for the disturbance. Unsteady critical layer 
analysis has also been applied to study the evolution of waves in stratified flows, Churilov 
& Shukman (1988), and also in the flow over compliant surfaces, Gajjar (1991). 

In this paper we consider the spatial/temporal evolution of an oblique mode propagat- 
ing in a two-dimensional compressible boundary layer with non-insulated wall conditions. 
In Part II of this paper Gajjar (1992) the extension of this work to the insulated case for 
0(1) and also small Mach numbers is also described. In fact the generic problem there is 
essentially the same as that obtained here and can be retrieved also by taking a limit of 
the current problem. Long waves, based on the Rayleigh scalings, are considered and the 
growth rates are assumed to be small so that the ideas of nonlinear critical layers can be 
used. The mode is taken to grow slowly in the direction of propagation of the wave, and 
the cross variations are taken to be small and negligible. The assumptions made imply 
that the analysis here applies only to the evolution of the first mode, in the terminology 
of Mack (1984). The higher modes have much larger phase speeds and stem from regions 
of locally supersonic flow which require a modification to the present analysis. 

In the following sections we analyse the disturbance quantities and obtain the gov- 
erning nonlinear unsteady critical layer equations describing the evolution of the mode. 
It is shown below that the evolution is described by a pair of coupled equations for the 
disturbance vorticity and density, with the density acting as a source of the vorticity. 
These equations are essentially the unsteady counterparts of those obtained in I. The nu- 
merical solution of these equations using a novel method based on Fourier decomposition 
in the spanwise direction and Chebychev collocation in the normal direction is described 
also. The same notation and nondimensionalisation of I is used and in particular (r,y, z) 
denote the streamwise, normal and spanwise non-dimensional coordinates, (u,v,w) the 
corresponding velocities, t is the time, (p, p, T), are the nondimensional pressure, density 
and temperature respectively. The Reynolds number R is taken to be large throughout, 
Moo is the free-stream Mach number, and F is the ratio of specific heats. The basic flow is 
taken to be a flat plate or pressure gradient boundary layer with wall heating and cooling, 
see I also. 


2. Problem Formulation 

It is convenient to introduce Squire coordinates (£, z) along and normal to the direction 
of propagation of the wave, with 




d_ 

dx 
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T,- Ri w4( + h M )+ah ^’ 

l^h’RH-oc-hr)^, 

and £ = R$h(ax + (3z - acht),z = R.\h(az - f3x), and X = Rhh 2 (ax + flz), with the 
wave growing on the slow X scale. The scale factor h,(k « 1) is introduced so that we 
are considering long waves on the Rayleigh scalings. The growth rates are thus of 0(h 2 ). 
Also c is the phase speed of the wave, <t is taken to be real and a — cos 8,(3 = sin 9, where 
9 is the direction of propagation of the wave. With respect to the Squire coordinates we 
have the corresponding velocities u = au + (3w,w = —(3u + aw. The basic flow is taken to 
a two-dimensional boundary layer flow with 

(u,v,w) = (ir B (x,y),o(R-l),0),p = PB = 1/rMl, 

T — T B (x y Y),p = ps(x y Y) = 1/Tb , 
where y = R~ 3 Y and we assume also that 

U B ,p B , —> 1 a* y-*oo, 

U B ~ XiY + \ 2 Y 2 + ... as y-»0+, 

Pb ~ flg 4* Rj Y + ... as Y — > 04*, 

Tb ~ 5o 4* S\Y 4* . . . as Y — * 0 4* • 

Small disturbances of size 6 are introduced and expansions for the disturbance quantities 
are considered in each of the regions, see Fig. 1. The disturbance size S = 0(h 3 ) has been 
anticipated below to allow the nonlinearity and unsteadiness to appear in the critical layer. 
Firstly in the region Y = 0(1) the various flow quantities expand as 

u = aU , b 4- 4- h* 4* • • • , 

v = h*v i 4* h 6 v 2 , 

w = — /3U b 4* h 3 wy 4- h*v>2 4* • •• , (2-1) 

P = PB 4- h*pi 4- h s p2 4- , 

P=PB 4- h 3 p! + h*p 2 + ..., 

and these are taken to be independent of z to this order. Substitution of (2.1) into the 
Navier- Stokes equations yields at leading order the solutions, 


4 



ill = AqocUby, vi = -aA 0 (U B , Pi = A, 

Pi = A 0 pbk> w = -PAqUby, 


(2.2) 


where A 0 = A(X)e* 7 ^ + c.c. , Pi = P(X)e‘ 7 ^ + c.c.and A, P are slowly varying amplitude 
functions, and 7 is the wavenumber of the oblique mode. At the next order the solutions 


t >2 


-"-il 


Y PH 


dy -f caAo^ + A\Ub — 

Y 


and 


PbUI 

Pi = P 2 - / aUspBVit dy, 

Jo 

oUbP if 


Pi(oYUb 

rps 


(2.3) 


«2£ = — ~ V 2 Y ~ 


Tps 


are obtained and A 2 , Pi are unknown function of £ and X. 

In region Z2 where K = h?, the properties of the basic flow in (2.1) and the solutions 
(2.2), (2.3) imply the expansions 


u = a(h\iY + h 2 \ 2 Y 2 ) + h 3 U! + h*u 2 + ..., 


v = h 5 v 1 + h a v 2 + . . . , 

w = — p(k\iY -f- h 2 X 2 Y 2 ) A h? w 1 -f- h 4 w 2 

P-PB + h 4 pi + h s p 2 + . . . , 

p = Ro + hRiY + h 2 R 2 Y 2 + h ? pi -J- h 4 p 2 -f- . . . . 
Substituting into the equations and solving gives 


v\ = —aAo^XiY, ui = aAiA 0 , 


Wi 


PX 2 YA 0 _ RiXiYAo 
(X iY-cY Pl (X 1 Y-c) 


<TU — 0 CCU\X 
V2 = aX, 


aAiA 2 Ao ( r _ 2 21 , 

— -W -c 2 } + 


PnRic 

arPpAj 


(P 2 t+Pix) occA 0 (. 2X 2 Pi 


P 0 a 0 Ai 


Ai xtg 


(2.4) 


where 77 = (AiF - c). In the above we have anticipated the condition that the normal 
velocity is zero at the wall which leads to relation (from the solution for vi) that 


These solutions show that the density, temperature and spanwise fluctuations become 
algebraically singular in the critical layer. The terms are jump conditions arising from 
the continuation across the critical layer at fj = 0. Before discussing the critical layer 
equations it is convenient to consider first the outer potential flow region Z3 and derive 
the amplitude equation describing the evolution of the mode. 

In region Z3 with Y = Y/h we have the expansions 

u = o + h 4 u i + h 5 u 2 + . . . , 
v = h 4 vy + h 5 uj -4- • . . , 
to = —/3 + h*w\ -f- h 5 W 2 + • ■ • , 
p = 1 h> py "l" • « • , 
p = Pb + h 4 pi + h*p2 + ■ • ■ • 

When substituted into the Navier-Stokes equations the equation for pi is just the Prandtl- 
Glauert equation which yields 

Pi = Pme-^e^, = ( 2 - 5 ) 

where reed quantities are assumed, Pi o is an unknown function of X and 

fl = (1 - r?-)* = (1 - COS 2 OMl,)*. 

IpB 

At the next order the e* 7 * component of pa , say Pai is given by 

p 31 = + YPyoxi^e-^ + p2ie~^ , 

Ad V 

and the e*^ componenet of tJj, Vyy satisfies 

icnVn \t„ = CT(1 2 n n ’ > ^ +70ftl- (2-6) 

3. Derivation of the amplitude equation 

The amplitude equation is derived by matching between the different regions and also 
imposing the condition that the normal velocity is zero. At leading order from (2.4) and 
(2.5) we get 

Py 0 = a 2 XyAR 0 c=^^. (3.1) 
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which gives the dispersion relation 


7 = ClXiRoc = Aii2oc(l — cos 2 . 


(3.2) 


Denoting the e* 7 ^ components of p 2 >-4i, by P 21 ,Au and A<j>f respectively matching 
of the normal velocities and pressure between Z1 and Z3 gives using (2.3) the relations 


where 


and 


V 2 i|y =0 = —Pi oh + iyacA + A u , 

a 


P 21 = P 21 - a 2 7 2 Af 2 , 


T f °° r 1 1 , 1 ,Ri 2A 2 , „ , 

1 J 0 [ p B U 2 B RoXjy* + R 0 Xjy ( R 0 + A a } 1] dy ' 


(3.3) 


h= r{ PB Ul-l)dy. 

Jo 

Matching the normal velocities in Z1 and Z2 gives from (2.3-2.4), 

*tAoA 2 , j x _ iac7i4 / 2A 2 , Pi xx 1+ , 2iX 2 ajcA 
r An A 1 = ; (~r 1- -5 ")Ai <f>i + 


ocR 0 Af 


Ai \i Rq 


Ai 


(3.4) 


Next setting the second order component of the normal velocity in Z2 to zero at the wall 
gives from (2.4), 


iD\A — acAx 


RqOLqX 


*7 - iyac 2 , 2A 2 Ri . . _ - 

P\*x - TT —P 21 + — (-T- + -4 = 0. 


U0O1A] 


A a v Aj R 0 


(3.5) 


where 


yRiC P 10 yc 2 a 2X2 , Hk. . . 

fli =^ + ^p?<T )+ Tr ( Ar + ^ )ln|c| - 


Finally eliminating P 2 i>-4 n using (2.6), (3. 1-3.5) leads to the amplitude equation 
. „ 7 7 t7^4ac 2 .2A 2 Pi,.., , , 

iD*A - 2 caA x - -4 (-1 + -±)(<j>+ - <f>~) = 0, 

Al ItQ 

where the D{ are given by 

Cy(l + ) ( PjQ ^ 3r> 2r f 2 2 , 2r / ^ 10 \ 

= — 2 fT — (“X' ~ 7 J2 + 7 a c + 7 2 i 1 (-=~), 

n 1 r n 2a 2 7 2 A 2 c 7 2 A 2 ,Pi 0vi 

a| 


(3.6) 


7 


and 


D 4 = D 1 - 


7 

JZoaoAi 


D 3 . 


Note that D+ is a real quantity. 

The jump conditions {4>i are determined from the properties of the critical layer 

and this is considered next. 


4. The critical layer equations 

Expansions for the critical layer follow readily from the solutions and expansions in 
Z2. In the critical layer with Y = ^ + h 2 r], we have 

H = a[hc 4- 4 — )] + h?U\ h*U 2 4" • • • > 

v = h*V 0 + h*Vi + h 7 V 2 + ... , 

-tU = — fihc 4- h 2 Wi 4- . . . , 
p — pb 4- h* P\ + h 6 P 2 + . • • , 
p — R 0 4* h~r h + .... 

Ai 

Substitution into the equations shows that the pressure is the same as outside the critical 
layer. The leading order solution for the normal velocity yields 

UoaAiVo = — Pit* 


The first nontrivial solutions lead to the critical layer equations 


Lir 1 = 0, 


LQ 2 = 


*1 tfPu 

R l ’ 


where the Q = U 2v and the operator L is defined by 

. A2C 2 <r.d d d 

L = a(A„ + + V,^. 

In order to match with the solutions outside we need 

wi ~ R\Tf 4- Re(Ae 1 ' 1 ^) + . . . , 

*17) 

Q 2 ~ 2a\ 2 T, + — 4 - ^)Re(Ae*t) + . . . , 

rj Aj Ko 


(4.1) 


(4.2) 
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as r} — > ±oo. Here Re denotes the reed part. Hence the evolution of the mode is governed 
by the solution of the equations (4.1) together with the amplitude equation (3.6). Using 
(4.2) the amplitude equation can be written in the form 


• 2 /•— /• oo 

_ vfc r y r 

2 -k\\ J 0 Jcc 


iD^A — 2 cxcAx 


Q 2 e-^ dr)d£ = 0. 


(4.3) 


The D 4 term in (3.9) and (4.3) gives a correction to the wavenumber and can be in effect 
be removed by redefining and is thus set to zero. Note that this also fixes <r and hence 
the correction to the frequency. 

A more standard form of the equations can be obtained by renormalising with the 
scalings 

7 T! = Ridl H, X = - Xo, 

7^ = ^* -f £o> A = d^Ae *^°, Q 2 = d 2 Q, 

where 

c 3 ,2A 2 Ri. 3 .2c<fjxi 

d »-8Af ( 17 + ^ ) ’ 1_( ^ } ’ 

2acd5 /2A 2 , i2i x 

^ = + 


to give 


(£ + &)’ 

*+z*L-a#A*r)\ 


an „sn 


v dll 


Also the boundary conditions (4.2) and the amplitude equation (4.3) reduce to 

Q~(l-J)Z + + ... as Z -+ ±oo, 


„ „ i2e(Ae*«‘) 

II ^ Z “h - — — — as Z — ► i oo , 

Z 
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and 


(4.6) 


* - 2w *00 

Ajf = - - / Qe~*'di'dZ. 

* Jo J-oo 

The constants X 0 ,£o can be chosen to match upstream conditions and A -» fax 
upstream ( X — * — oo), where the scaled linear growth rate is now tt. 

The probl em (4.3-4.6) is similar to that obtained by Goldstein & Wundrow (1990) in 
their investigation of the stability of hypersonic boundary layers. Note that viscous effects 
have been neglected above which means that the scale factor h is such that h» 0{R-Z). 
If h = 0(R~ i«) (4.4) is modified to 


L(Q) = JRe(iAe « * )|| + + 7e 


d 2 n 

dz*' 


X( n )=7a^, (4.7ofe) 

where the operator L is defined by 

L = —=r + Z-^r Re(iAe l t 

dX df v ’ dZ 

and the parameters 7c ,7<j>7e are dependent on the viscosity law used. With the inclusion 
of viscosity the wavenumbers are in effect 0(R % ). We also have 


Tc = 


5 0 /i(5 , o)fi-U- 7 

7aAadJ 3 

7e=7dJ( 


7 d = -fc/Pr, 


Prf^{So) 

Rot*(S 0 ) 


- 1 ), 


(4.8) 


where Pr is the (constant) Prandtl number and is the viscosity law. The properties 
of the steady state or equilibrium version of (4.7) with j e = 0 are discussed by Gajjar & 
Cole (1989). 

It can be seen from the above scalings and (4.8) that as far as the oblique-wave is 
concerned, the 6 dependence can be scaled out of the invisdd problem, but not for the 
viscous problem. In fact the more oblique the wave the larger the effective viscosity for 
fixed values of the other parameters. 

For an insulated flat plate boundary layer the same normalized problem (4.4) is ob- 
tained, although the J parameter is the limit of that given here as R x tends to zero, see 
Gajjar (1992). We note also that Wundrow (1989) has obtained equations similar to (4.7) 
in his study of the stability of compressible mixing layers and he presents several numerical 
results. 
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Linear growth rate. 

For the linear problem the standard jump conditions of in across the critical layer 
when substituted into (3.6) determines the linear growth rate as 


K r 


'yen ,2X 2 f2i 
2Ai ' A, + R 0 } 


(4.9) 


A more useful form for the growth rate, wavenumber and frequency can be obtained 
by first recasting the various expression in unsealed form i.e. without the h factors present. 
From the multiple scaling used it can be seen that the effective wavenumber, frequency, 
and growth rate are given by 

wavenumber: R* 7* = (jh + ...), frequency: R j«, = R* (7 h 2 c cos 0 + . . .), 

growth rate: Ra k* = R* (« r h 2 -f- . . .), phase-speed: c, = — — — . 

7, cos 6 

For the flat plate boundary layer flow, taking the Prandtl number to be unity, the 
basic flow is given by in Dorodnitsyn-Howarth variables, 

U B (x,Y) = f ( ;,), (2x)^= f Y PB dY 

Jo 

where f(rj) satisfies the equation 


(fi(T B )PBf"y + sf = 0 , m = ao) = o,/'(oo) = 1, 


and 


— = Tb = (1 - h)U B +h + l( r - i)M^(t b + u B )( 1 - U B ). 

PB * 


Here tb is related to the wall temperature T\y by 

T«r = (1 + iwi(r - l))< t 

and fi(T B ) is the viscosity law. In particular tb = 1 corresponds to the case for insulated 
walls. 

The above asymptotic expressions for the wavenumber (3.2) and growth rate (4.9) 
thus reduce to 

(4.10) 


2 U>„Ro/"( 0 ), 2 2 mi 

7* = J V (1 ~ cos 2 6) > , 


u>,n ( 1 

K * = a l 

cos 6 \ 


V^cos 6 

-tb' 


tb 


d(A J* )pB ) Tw 

1 + «I& 

2p(T B )p B 


) 


(4.11a) 


,=o 


The expressions (4.10), (4.11a) give the leading terms in the asymptotic expansion for the 
wavenumber and growth rate in terms of the frequency in this long wavelength approxi- 
mation. It is seen that for heated walls, with <5 > 1 the growth rate is positive and the 
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growth rate increases with increasing wall heating. Also the above expression shows that 
for fixed frequencies the oblique waves have larger growth rates than the 2D modes. 

Clearly when = 1 i.e for the insulated case, the growth rate reduces to 0(h 3 ). The 
insulated case is considered in Part II of this paper Gajjar (1992), where it is shown that 
for 0(1) Mach numbers, the growth rate term is given by 


.(r-i)AC / 

Tw cos 2 07 * y 2 h(T b )pb /^_ 0 


7TW? 


(4.116) 


Further when the Mach number is much smaller of 0{h a ) then the growth rate reduces to 
0(h 4 ) and the A 4 contribution comes into play. In this case K m above is modified to 


7TU>: 


«. = 


COS' 


*7. 


(r - 1)M£ 

Tw 


(■* 


)- 


C* 


dit(TBlPBl Tw 

2 KTb)pb I 0) 3 


(4.11c) 


J TJ = 0 


The above expressions taken together show that for long waves, for insulated or isothermal 
walls, a generalised expression for the growth rate asymptotically valid to much higher 
order is 

*•= 2null^ BDUB) < 412 ‘> 

where D = d/dY and pb,T)Ub, D(pbDUb) are evaluated at Y = Y e . The expression 
(4.12a) reduces to (4.11) by expanding the various quantities in (4.12a) and noting that 
for small c*, Y c = c*/Ai to a good approximation. For neutral modes (4.11c) agrees with 
that given in I. 

From (3.2) a leading order expression for the frequency is given by 


= cl cosODUbPbO- ~ cos2 MC)$. (4.126) 


The growth rates as computed from (4.11), (4.12a) are compared with numerically 
computed values from the solution of the linear compressible Rayleigh problem* and the 

r 3 / 2 

results are shown in Figures 2-4. The Sutherland viscosity law p(T) = j— — — together 

(0.5 + T ) 

with T = 1.4 was used in these comparisons. The growth rates are plotted against the 
phase-speed because, for fixed frequencies the phase-speed as computed from (4.12b) would 
not be particularly accurate as (4.12b) is derived from just the leading order expression 
for the wavenumber. Using the phase-speed as a parameter and computing the other 
quantities from the various expressions written in terms of c» gives a much more meaningful 
comparison. It is seen, Fig. 2(a-d), that for the non-insulated case whereas the formal 


* The author would like to thank Dr. L. Hultgren for providing the code used in the 
computation of the inviscid spatial eigenvalues. 
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leading order asymptote (4.11a) works well only for low Mach numbers and very low 
frequencies, the generalised expression give excellent agreement for a range of frequencies 
and oblique angles upto moderately supersonic Mach numbers. Figures 3(a-e) show the 
comparisons for the insulated case also. The expression (4.11c) and the generalised form 
(4.12a) again give excellent agreement with the full solutions for the same range of Mach 
numbers. For increasing Mach numbers the critical layer moves away from the wall and 
the assumptions used in deriving (4.11), (4.12a) become less valid. Comparisons of the 
asymptotic expression for the wavenumber against numerically computed values in Figures 
4(a-f ) show that the leading order asymptote works well only for small c, and further terms 
would need to be incorporated for a better comparison. 

5. Numerical solution. 

The method used to solve (4.5-4. 7) is different from that used previously and is de- 
scribed below. The method was developed so that the calculations could be carried on 
the AMT DAP-510 machine which is a massively parallel SIMD system comprising 1024 
processors arranged in a 32 by 32 matrix. Since the individual processors of this machine 
are relatively slow, the machine is most efficient only when the underlying algorithm is 
highly parallel a feature of the current method. 

The vorticity Q and density II were decomposed into individual Fourier modes by 
writing 

<j = (i-j)z+ (if; <?.«*»<• +*.c), 

2 V 

n = Z + (t^ n „e i “f +c.c) 

0 

to give a set of equations for the Q n , II n similar to those given in Goldstein & Hultgren 
(1988). The system was truncated to solve for (N + 1) Fourier modes and a predictor 
corrector scheme was used to advance in time. The resulting equations at a particular 
time level (j + 1) are then of the form 

<?i +1 + ikZQi +1 - 7e<?izz = r„, (5.1) 

for each Fourier mode Qk , (Jfe = 0, ..., N) and the nonlinear contributions are contained 
in the r* term on the right hand side. These were then solved using Chebychev collocation, 
by first mapping the infinite interval into [-1,1] with the transformation 

Z = 2y/y/l-y* t 
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and writing 


M 

<?i +I = !>!"%>. 

J=0 

to give from (5.1) a set of linear equations for the Chebychev coefficeints a\ k \ The inverses 
of the coefficent matrices were pre-computed and stored, and the were obtained by 
doing matrix- vector multiplications for each Fourier mode. The method is thus fully vec- 
torizable/parallelizable and suitable for massively parallel systems like the AMT DAP-510. 

Note that since the infinite region is mapped into the finite region [-1,1] the Haynes 
(1985) procedure does not need to be used. Another advantage of the method is that full 
spectral accuracy in and Z is obtained. As described above the method requires the 
pre-computation and storage of (i\T+l) complex (M+l) by (M+l) matrices, since as (5.1) 
shows the differential operator involves an explicit dependence on k with the ikZ term. A 
modification of the method was used to avoid the large storage problems by writing (5.1) 
in the form 

Qi +1 + iNZQi +1 - JcQizz = i(N - k)ZQk + r k . 

This worked quite well in practice. 

One final point is the evaluation of the integral in the scaled amplitude equation 

(5.2) 

With the mapping the integral becomes 


Ajf = —i f Qi dZ. 

J — OO 



and was evaluated using Gaussian quadrature with the Chebychev poylnomials as basis 
functions. The quadrature points qj — cos(^^^), j = l,...,n were used in the eval- 
uation of the integral as this avoids difficulties at the endpoints y = ±1. This too can be 
implemented quite efficiently as a matrix vector multiplication. 

6. Results and discussion 

The numerical method was first tested by computing some test cases of the Goldstein 
& Hultgren (1988) problem for shear layers. ( This corresponds to taking H = 0, J = 0 in 
(4.7) together with an additional term on the right hand side of (4.7a), see equation 
(3.16) of Goldstein & Hultgren (1988)). Some sample results are shown in Figures 5(a- 
d) and these are in excellent agreement with those of Goldstein & Hultgren (1988) and 
Goldstein & Leib (1988). It was found necessary to use a large number of Chebychev modes 
in the normal direction, typically 150 modes, together with very small timesteps, especially 
for the inviscid cases, in order to reproduce the results. One reason for this is that as the 
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disturbance evolves downstream, the vorticity develops into the familiar cat’s eyes pattern 
and therefore many modes are required to resolve the very thin shear layers which form. 
The calculations were quite sensitive to the accurate computation of the integral in (5.2) 
and here 120 modes were needed for the later stages of the evolution of the disturbance. 

Some results for the current problem (4.7) axe shown in Figures 6-7. To reduce the 
number of parameters in (4.7) the Prandtl number was taken to be unity and the Chapman 
viscosity law was used giving -y e = 0. For J zero or small and positive, Figures 6(a-d), the 
results are similar to those of the shear layer calculations, Figure 5, see also Goldstein & 
Hultgren (1988). The growth rate of the wave follows the linear value closely and then 
decreases as the nonlinear terms become more dominant. The J — 0 case is also a special 
case of that arising in the instability of the flow over compliant surfaces, Gajjar(1991). In 
Figure 6d we show the effect of inputting incorrect initial conditions. The growth rate after 
rapid oscillations settles down to its correct linear value. For larger values of J positive, 
Figure 7 the growth rate decreases initially from the linear value but far downstream this 
is reversed and very large fluctuations are present. The effect of viscosity can be seen in 
Figures 7 which show that the smaller the viscosity parameter y c the greater the effect on 

the amplitude in the nonlinear regime. In figure with j c = 1 the disturbance reaches a 
large amplitude and then settles down with small oscillations about this amplitude. For a 
smaller value of -y c = 0.5, in Figure 7(c,d) it is seen that after reaching a peak amplitude 
the effect of nonlinearity and the forcing from the density causes very large fluctuations 
in the amplitude. Near X = 3 the amplitude drops rapidly and then increases again with 
the decrease/increase taking place within a short X-scale. 

The computations in Figures 6,7 do not go far enough in X to be able to ascertain 
what the the limiting properties of the downstream evolution axe, although judging from 
the work of Goldstein & Hultgren (1989), it is expected that the properties of equilibrium 
compressible critical layers as dicusssed in Gajjar & Cole (1989) do come into play at 
some stage with outer diffusion layers also present. This however requires further analysis 
and more extensive numerical computations of the nonlinear problem. The nonlinear 
calculations require considerable resources, especially the computations with small viscosity 
where very small timesteps are necessary to allow a stable calculation. 
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Figure 2d 

Figures 2(a-d) Comparisons of the imaginary part of the waveneumber 7 * versus the 
phase speed c* from the full numerical solutions of the linear Rayleigh problem (solid 
line), with asymptotic results, see the text. The thicker dashed line is the generalised 
expression (4.12a), the thinner dashed line the asymptote (4.11a) for the non-insulated 
case. The parameter tb = 1 . 1 . M = Moo is the Mach number and 6 the oblique angle. 
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Figure 4f 

Figures 4(a-f) Comparisons of the real part of the wavenumber from full numeri- 
cal solutions of the Rayleigh problem (solid line) with the asymptotes (4.12b) (thicker 
dashed line) and (4.10) (thinner dash). Other parameter as in Figure 2. 








Figure 5(a-b) (a) A plot of the scaled wave amplitude A and (b) the scaled growth 

rate Re(A' /A), as a function of the the spatial coordinate X for the Goldstein & Hult- 
gren (1988) problem (see text) with = 1 and U = 1. 

Figure 5(c-d) As in Figure 5(a) except 7 C = 0, and U = 3. In Fig. 5(d) the real and 
imaginary parts of — A 1 /A tire plotted. 
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Figure 6d 

Figure 6(a-b) Solution of the nonlinear problem (4.7) and the figures show the evo- 
lution of the scaled wave amplitude as a function of X. The different curves are for j c = 
10, 1 and 0.5 with J = 0. 

Figure 0(c-d) As in Figures 6(a-b) except J = 0.1, 7 C = 1. 


28 







J-0.5 7 e * 1.0 



-2 0 2 46 

X 

Figure 7b 


29 








REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of Information Is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of Information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

March 1993 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Nonlinear Evolution of the First Mode Supersonic Oblique Waves in 
Compressible Boundary Layers. Part I-Heated/Cooled Walls 


6. AUTHOR(S) 

J.S.B. Gajjar 


WU- 505-62-21 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-7718 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 106087 
ICOMP-93-08 



11. SUPPLEMENTARY NOTES — 

J.S.B. Gajjar, Mathematics Department, The University of Manchester, Manchester, M13 9PL, England, and Institute 
for Computational Mechanics in Propulsion, Lewis Research Center (work funded under NASA Cooperative Agree- 
ment NCC3-233). ICOMP Program Director, Louis A. Povinelli, (216) 433-5818. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 34 


1 3. ABSTRACT (Maximum 200 words) 

The nonlinear stability of an oblique mode propagating in a two-dimensional compressible boundary layer is consid- 
ered under the long wave-length approximation. The growth rate of the wave is assumed to be small so that the ideas 
of unsteady nonlinear critical layers can be applied. It is shown that the spatial/temporal evolution of the mode is 
governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity and density. Expressions for 
the linear growth rate show clearly the effects of wall heating and cooling, and in particular how heating destabilizes 
the boundary layer for these long wavelength inviscid modes at 0(1) Mach numbers. A generalized expression for 
the linear growth rate is obtained and is shown to compare very well for a range of frequencies and wave-angles at 
moderate Mach numbers with full numerical solutions of the linear stability problem. The numerical solution of the 
nonlinear unsteady critical layer problem using a novel method based on Fourier decomposition and Chebychev 
collocation is discussed and some results are presented. 


14. SUBJECT TERMS 


Critical layers; Compressible; Boundary-layer stability 


15. NUMBER OF PAGES 

32 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 18. S EC URITY CLASSIFICATION 19. SECURITY CLASSIFICATION 

OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


























